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ABSTRACT 

We show that the evolution of the quintessence energy density CIq is model independent in an accelerating universe. 
The accelerating behaviour has lasted at most 0.5 e-folds of expansion assuming a present day value of fig = 2/3. 
The generic evolution differs from the exact solution by a percentage given by r = 27q, with jq = (pq+Pq)/ Pq- For 
a small jq the evolution remains a good approximation for a long period. Nucleosynthesis bounds on J7q suggest 
that the model independent solution is valid for at least 12 e-folds of expansion which includes then the scale 
of radiation and matter equality. We can, therefore, establish model independent conditions on the cosmological 
parameters at some factor scale Oj <C a . Finally, we discuss the relevance of this result to the cosmological 
coincidence problem. 
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Recent observations show that the universe has entered an accelerating expansion regime Q. If these 
observations are confirmed it would be the first experimental evidence for an energy density different from 
the radiation or matter and with negative pressure, i.e. for a " cosmological constant". The cosmological 
constant can be given by constant energy density with equation of state p\ = —p\, where p\ is the 
pressure and p\ the energy density, or it could be parameterised in terms of a slow varying scalar field 
Q, quintessence, with a potential V(Q) and a time varying equation of state pq = (jq — 1)pq, with 
PQ = \Q 2 — V(Q) and pq = \Q l + V(Q). The equation of state has < jq < 2, it is model dependent 
and it may vary with time. For an accelerating universe we have the condition (jq — 1)Oq < —1/3, 
where £1q is the ratio between the energy density of Q and the critical energy density. 

An important question is why the energy density of the cosmological constant is of the same order of 
magnitude as matter even though the expansion rate of both energies is quite different today. This is 
the coincidence problem Q. A possible answer is by setting the initial conditions such that the Q field 
start to dominate only recently. However, this solution introduces a fine tuning problem on the initial 
conditions. To avoid such a fine tuning problem tracker solutions have been proposed Q. In this models 
the energy density of the scalar field Q redshifts mimics the dominant energy density (radiation) at the 
beginning with 7q > 1 and when matter dominates the equation of state for Q changes to jq < 1. With 
this kind of fields one avoids the initial fine tuning problem but they have difficulties to explain (without 
any fine tuning of the potential) the central values of £7q = 2/3 db 0.05, 7q < 0.35 ± 0.07 || since they 
have 7q > 0.3 J2|. However, they remain a very interesting possibility. Other models with more the one 
term in the potential V = V\ + V2 where one term dominates during radiation and scales as radiation 
or matter while the other dominates at a recent time giving the acceleration of the universe have been 



studied [12]. In particular if during the radiation dominated epoch on has a potential V\ = aQ 4 , it will 
redshift as radiation and the ratio p r / pq remains constant. On the other for V\ = aQ 2 the ratio p m j pq is 
constant. In either case the initial value can be set to p r /pQ = 0(1) or p r /pQ = 0(1) avoiding an initial 
fine tuning problem but not the coincidence problem. Finally, we can invoke the anthropic principle to 
try to explain the coincidence problem where the likelihood of our present day cosmological parameters 
are determined 

Different models have been proposed to give an accelerating universe and a general analysis can be found 
in ||. We have the exponential potential, V ~ e _A< 2 with A < 3 H,|l]|, the inverse power potentials 
V ~ l/(J n , n > V ~ e^ Q II or V ~ (c + (Q - Q ) n ) a e- XQ Q and mixtures of any of these 

potentials. All these potentials have a finite A = — Vq/V < 0(1), with Vq = dV/dQ, when V approaches 
its minimum and this is indeed required for any model to give an accelerating universe since in this 
limit the potential energy dominates the energy density of Q. Some of these potentials may arise from 



non-perturbative effect or form string theory [14], Jig ]. 

In this work we address the reverse coincidence problem. This is, we start by imposing present days 
values of Oq, £l m and and we go backwards in time. We will show that the evolution of Q,q while the 
universe accelerates is model independent. We can therefore establish generic initial conditions for Qq at 
the beginning of the accelerating epoch. Furthermore, as long as 7q is small (say 7q < 0.4) the evolution 
of £1q remains model independent (within an error of r = 27Q = 20%) and this allows to set the initial 
conditions at a much smaller scale factor a% compared to today's value a a , e.g. at radiation and matter 
equality ai/a ~ 10~ 3 or at nucleosynthesis a,i/a ~ 10~ 10 . The range of validity is model and initial 
conditions dependent, however, in most models the universe have a scaling regime with 0,q <C 1, 7q <C 1 
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before the acceleration epoch leading to a large number of e-folds of model independent expansion. This 
regime seems to be necessary due to the condition Qq < 0.1 at nucleosynthesis @. Here, we do not solve 
the coincidence problem but we contribute to establish model independent behaviour of the universe at 
the accelerating era and some e-folds before that, say a m _i_. The last step would be to mach those initial 
conditions at a m .j. with a particular model coming from earlier times. 

In a spatially flat Friedmann-Robertson- Walker (FRW) Universe, the Hubble parameter is given by 
H 2 = \p = \{pb + Pq) where the barotropic fluid is described by an energy density pj, and a pression pb 
with a standard equation of state pb = (76 — l)p&, with 7;, = 1 for matter and 75 = 4/3 for radiation. 
Since we will study the evolution starting from today, we will take 75 = 0, however we will leave the 7& 
dependence in the equations. The interaction between the barotropic fluid and Q is gravitational only 
(we take 8irG = 1). In terms of the of the variables x = Q/^/6H, y = y/V/3H [l^] the cosmological 
evolution of the universe is given by 



x N 



VN 



-3x + ^Xy 2 + ^x[2x 2 + 7j,(1 - x 2 - y 2 )] 

-\J\xxy + \y\2x 2 + 76 (1 - x 2 - y 2 )] (1) 
H N = -l H [ lb {l-x 2 -y 2 ) + 2x 2 \ 

where N is the logarithm of the scale factor a, N = ln(a), Jn = df/dN for / = x,y,H and X(N) = 
— Vq/V. Notice that all model dependence in eqs.(|l|) is through the quantities X(N) and the constant 
parameter 7^ and generic solutions can be obtained ||,||. The flatness condition gives 1 = Qq + fi& and 
from eqs.(^) the evolution of Qq = x 2 + y 2 and 7q = 2x 2 /(x 2 + y 2 ) is 

(n Q ) N = 3(76 -7<j)n<j(l - Oq) 

(7q)jv = 7q(2-7q)(-3 + /|^ q ) (2) 

An accelerating expansion universe requires p + 3p = Pb(3jb — 2) + Pq{3^q — 2) < which gives the 
constraint 

( 7Q - 1)0 Q < -1/3 (3) 

where we have taken the barotropic fluid as matter, i.e. 7^ = 1. Equivalently, we can write condition (|3[) 
in terms of x as x 2 < (Qq — 1/3) /2 and at present day this sets an upper value of x 2 < 1/6 and 7q < 1/2 
for today's central value of 0, o q = 2/3 (from now on the subscript o refers to present day quantities). If 
we take the preferred value 7 g q < 0.35 the condition on x is x 2 = \£l Q loQ < 0.11, which is slightly 
smaller then the previously obtained one, and y 2 = £1 q(1 — 7 q/2) > 0.55, i.e. we have y 2 > 5x 2 . The 
solution to eq.(|2|) is only model dependent through the quantity X(N). Solving eq. (Q) for J7q we obtain 

n Q (N) = ^ = (4) 

i-n oQ + n oQ e-^ N °- N ^ 3 f^ dN 

where 0,q(N o ) = Q q q and we have N < Nq for a time t < to. In the limit of S7 q -C 1 and 7q ~ cte eq.(H) 
reduces to fig ~ e 3N ^ b ~ 1Q \ However we are not interested in this region since the initial condition has 
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VLoQ ~ 2/3. During the accelerating period we have x 2 < y 2 and 7g < 7^ and we can expand eq.(f|) as a 
function of R = 3 / 7gdiV. Eq.(Q) becomes 

"^M i+ i-^X^-"> r+o(r2) ) <5) 

with 

nQ = l-n oQ + Q oQ e~^o~N)- ( 6 ) 

Notice that as long as R is constant or R <C 1, eq.(|5|) is model independent and it is certainly valid during 
the acceleration period of the universe. We see, therefore, that the accelerating expansion of the universe 
is model independent and the energy density of quintessence evolves as £Iq. Furthermore, fig remains 
a good approximation to eq.(Q) as long as 7q is small. We can establish a relationship between the 
amount of discrepancy that we wish to accept between Qg and f2g. For a r percentage of difference, i.e. 
Oq/Oq = 1+r we have r = — — — r ~ ^> ( at i arge jv _ j\r). Since {^q)n is proportional to 

7q it varies rapidly when 7g is not small as can be seen from eq.(||). R becomes large when 7qtv — ~6jQ 
(as a first approximation) and we find R = 3j "fQdN = (A , jq)/2. If we allow an r = 20% discrepancy 
between the exact and the model independent fig, this will happen at A7Q = 2r = 0.4. Since before 
the accelerating epoch one has x <C y and 7g <C 1, the model independent energy density stops being a 
good approximation only until jq ~ 2r. When will this happen is model and initial condition depended, 
however, for a large number of cases Nq — Ni can be quite large (iVj is the beginning of the model 
independent evolution of Q). In fact, if ,for example, we set initial conditions such that x < 2y, i.e. 
7q < 0.4, at Ni (regardless of its value) the evolution of f2g is model independent (within a maximum 
of 20% discrepancy). 

In fig. la and lb we show the behaviour of the exact f2g by solving the dynamically eqs.(||) numerically 
compared to f2g for the same exponential potential V = Voe -3 ^/ 4 but with different initial conditions 
x = 0, y = 0.001 and x = 0.5, y = 0.001 respectively at N = 7^ N a . To appreciate more the difference 
between these two quantities we plot Qq/VLq and we include in the graph the acceleration parameter 
a = p + 3p/(pf, + 3pb) H and 7q. We observe from the graphs that the concordance between Qq and 
Qq is remarkable as long as the universe expands in an accelerating way (a < 0) and 7q remains small. 
We would like to mention that running the models down from present day the solution is very sensitive 
to the initial conditions at N . This is because we have a late time attractor and models with small 
variations at a late time may have come from quite different values at AiV large. We also show in fig.lc 
two different models (V — 1/Q and V — e -1 /^) with initial conditions x = 0, y = 0.001 at N = 0. 

Form fig.l we see that as long as x, 7q are small the agreement between £Iq and is very good. In order 
to avoid a rapid increase in Qq (going backwards in time), which would contradict the nucleosynthesis 
bound (Qq < 0.1), we require precisely this kind of behaviour. So, we expect our £Iq to be a good 
approximation for a long period of expansion that could go as far as matter and radiation equality or 
nucleosynthesis. Furthermore, if we run the models backwards the "early" time attractor has fig = 1 
(with x = 1, y = 0) and the scale at which the transition between 7g > 1 and 7q < 1 takes place is a 
symmetric point in Q,q as a function of N. Therefore if Qq <C 1 for AiV = 2M then our generic solution 
will be valid for at least M e-folds of expansion. If we assume that the quintessence energy density has 
been increasing since nucleosynthesis (i.e. no extra complications) then fig is valid for at least 12 e-folds, 
i.e. it includes matter and radiation equality scale. 
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Figure 1: In fig.la and lb we show the behaviour of Qq (solid line), 7q (dash-dot), a — p + 3p/(pb + 3pb) (short dashed 
line) and £Iq/Q.q (long dashed line) for an exponential potential with A = 3/4 and for two different initial conditions at 
N = 0, x = 0, y — 0.001 and x = 0.5, y = 0.001 for graph la and lb respectively. In fig.lc we show the behaviour of £Iq, 
7q and the ratio Qq/CIq for a potential V ~ e~ x ^ (thin solid line, short dashed and long dashed lines respectively) and 
for a potential V ~ 1/Q (thick solid line, dashed-dot-dot and dashed-dot line respectively) with the same initial conditions 
x = 0, y = 0.001 at N = 0. Notice that in all cases the model independent solution agrees quite well with the exact solution, 
for a large number of e-folds and the discrepancy becomes significant when 7q starts to increase. Finally, in fig. Id we plot 
the Q,q vs N. 
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We can invert eq.(||) to get the number e-folds of expansion AiV = N a — JVj (see fig. 2b) as a function of 

n iQ (Ni)/n oQ (N ), 



AN = -Log 
3 



QoQ / 1 — fii 



Log[l + R] 



(7) 



We can therefore determine the amount of quintessence On at different earlier times. For example, the 
number of e-folds of accelerating expansion is at most AiV = 0.46 since Oq must be greater or equal 
than 1/3 (cf. eq.@). During structure formation with <2j/a ~ 1/2.2 we have 0,q = 0.15, at matter and 
radiation equality (useful for tracker models) ai/a Q ~ 10 -3 we get Q,q = 2 x 10~ 9 and at nucleosynthesis 
with ai/a ~ 10~ 10 we have = 2 x 10~ 23 . With these results we can establish the initial conditions 
from which the model independent expansion is valid and we only need to concentrate on the behaviour 
at earlier times. 

To conclude, we have seen that the expansion of the accelerating universe is model independent and it 
is given by £Iq. As long as jq is small the evolution of £Iq remains a good approximation to the exact 
solution and the percentage difference is given by r = 2"fQ- Due to the nucleosynthesis bound (Qq < 0.1), 
we do not expect (going backwards in time from today's value) CIq to have a decrease and then a rapid 
increase and therefore our solution should be valid for at least 12 e-folds of expansion. We can then 
impose model independent condition at much smaller scales then our present day one, like at radiation 
and matter equality. This result contributes to solve the reverse coincidence problem since we can trace 
back the problem of initial conditions to a scale in which the universe expanded in a non-accelerating 
way and possibly as matter or radiation. 

A.M. research was supported in part by CONACYT project 32415-E and by DGAPA, UNAM, project 
IN-103997. 
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